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1 Introduction 



The aim of this article is to describe a new method to compute exact representations of time- 
dependent correlation functions in quantum integrable lattice models. For that purpose, we 
consider the example of the XX Z spin-^ Heisenberg chain in an external magnetic field. For 
simplicity, we will here mainly focus on the dynamical correlation function of the third compo- 
nent of spin. It will be clear however that the method is general and can be applied to other 
cases as well. 

This work is a continuation of PP, where we have derived a master equation for the (time- 
independent) correlation functions of the XXZ chain. In the end of that article, we have 
announced a multiple integral representation for the dynamical correlation function. We 
give here a proof of that result. 

The Hamiltonian of the periodical XXZ spin-^ Heisenberg chain in a magnetic field |2j is 
given by 

H = H^^^-hS^, (1.1) 

where 

M 

jj(o) ^ ^ (cr^fT^^+i + CFl,crl^+i + A(cr^C7^_^i - 1)) , 

m=l 
1 ^ 
m=l 

Here A is the anisotropy parameter, h denotes the external classical magnetic field, and a^'^ 
are the spin operators (in the spin-^ representation) associated with each site of the chain. The 
length of the chain M is chosen to be even. The simultaneous reversal of all spins is equivalent 
to a change of sign of the magnetic field, therefore it is enough to consider the case /i > 0. The 
quantum space of states TL is H = ®^^iTLm: where Tim ~ is called local quantum space. 
The operators Cm^'^ act as the corresponding Pauli matrices in the space Tim and as the identity 
operator elsewhere. 

The time-dependent local spin operators are defined as 

fT^y'^(t) = e^-^V^^'^e"*^*. (1.2) 

Since [cr^, Sz] = 0, we have for the local operator of the third components of spin 

(^m[t) = e CTmf^ ■ (1-3) 

Hence, the dynamical two-point o"^ correlation function at zero temperature is given as the 
following mean- value: 

9zzim,t) = (V^g|afe^^'°'*<+ie-^""*|^3), (1.4) 
where denotes the ground state of the Hamiltonian Hl.lj) . 
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The method to compute eigenstates and energy levels (Bethe ansatz) of the Hamiltonian 
was proposed by Bethe in 1931 in [2] and developed later in I^Elini- The algebraic version 
of the Bethe ansatz was created in the framework of the Quantum Inverse Scattering Method by 
L.D. Faddeev and his school [TllHlin]- Different ways to study the time-independent correlation 
functions of this model were proposed in the series of works (see e.g. ^1 ^1 El El El 
El El El El Uni 1^ ) ■ As for the dynamical correlation functions, up to now, the only known 
exact results concern the case of free fermions A = 0|lllEllllEllinillZlIlHl[ini- 

Our aim is to obtain a multiple integral representation of the expectation value ()1.4|) for 
arbitrary A > — 1. For this purpose we use two different approaches leading (hopefully) to the 
same answer. Both of them are based on the algebraic Bethe ansatz. 

The first method ^1 El Ej consists in acting with the operators af and c^+i (*) on the 
ground state {^pg\ to produce a new state, say {ip{z,m,t)\. Then one can compute the resulting 
scalar product {tp{z,m,t)\tpg). In the time-independent case this can be done by the algebraic 
Bethe ansatz and the explicit solution of the quantum inverse scattering problem for the o"^ 
operators in site 1 and m + 1 (see jl3l l3Uj). However, for the calculation of the dynamical 
correlation functions, we also need to compute the action of e*-'^'"'* on arbitrary states of the 
chain. To achieve this we use the fact that the Hamiltonian H^^^ can be constructed as the 
logarithmic derivative of the quantum transfer matrix T(A) : 



H(0) = 2sinh«^!I^r-i 



2Mcosh?7, with cosh?7 = A. (1-5) 



^-2 



Hence, from the Trotter type formula |3H I32j. we can obtain the evolution operator e^^^°^^ as 
the following limit: 

±..(H(o)+2Mcosh,) ^ hm (T(^- + e]-T-U^-]]^\ e = l 



g±^n/i-+^McosM] = ^im^ (^(i • '^^ iijj ' ^ = f 2itsinhr/. (1.6) 

It is moreover possible to express T^^{r]/2) in terms of T(— ry/2). Therefore the action of e*-'^^"'* 
reduces to the one of some product of transfer matrices. This enables us to apply directly the 
results of JSj and where we computed the action of any product of transfer matrices on 
arbitrary states in a compact form. In the present case it leads to a master equation for the 
dynamical spin-spin correlation function gzzi'm,t), or more precisely for its generating function 
(see Theorem 13.1(1 . 

The second method to handle Qzzi^, t) is to insert between the two operators a sum over 
a complete set of eigenstates of the Hamiltonian, leading to 

gzz{m,t) = Y.{^M\i){iWUMe'^'''-'''^'^ (1-7) 

i 

where = Ei\i) and H^'^^l^jg) = Eo\i)g). In P we have shown that for the time- 

independent case the master equation method allows one to define such a sum in a precise 
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manner. It uses in fact a twisted version of the transfer matrix T(A) and the knowledge of 
determinant representations for the form factors (e.g., corresponding to the matrix elements 
{'^gWi\i))- This method can easily be generalized to the time-dependent case. So, using the 
technique developed in we are able to re-sum completely the form factor expansion for the 
dynamical correlation function and to show that it leads indeed to the time-dependent master 
equation obtained by the first method described above. 

This time-dependent master equation gives a representation of the dynamical cr^ correlation 
function of a finite XX Z chain. In the thermodynamic limit M oo, the model exhibits three 
different regimes depending on the value of the anisotropy parameter A and on the magnetic 
field h. In this limit, following |13 | I14 | [T5 | IT], we obtain a multiple integral representation for 
the dynamical correlation function in the massive (A > 1) and massless (—1 < A < 1) 
regimes. It is given as a sum of multiple integrals (see (|f).16|) ) involving the time dependence 
through the factors e*^*, where E is the bare one-particle energy. 

This paper is organized as follows. In Section 2 we introduce some useful definitions and 
notations and recall briefly the properties of the twisted transfer matrix. The time-dependent 
master equation (|3.8() for the generating function of the dynamical correlation function is 
given in Section 3. In the following two sections we present two independent proofs of 1)3. 8() . 
The first one uses the Trotter-Suzuki procedure for the explicit calculation of the action of 
^iH<-°h arbitrary state of the chain (Section 4). In Section 5 we derive the time-dependent 
master equation via the form factor type expansion. In Section 6 we obtain a multiple integral 
representation for the dynamical correlation function in the thermodynamic limit. 

Some technical details are gathered in two appendices. In Appendix ^ we prove several 
lemmas on the properties of the solutions of the twisted Bethe equations and of the correspond- 
ing eigenstates of the twisted transfer matrix. These lemmas are used in Sections 4 and 5. In 
Appendix El we compare the representation 1)6. 16() for the dynamical correlation function at 
A = with the one obtained in ,29' in the case of free fermions. 

2 Algebraic Bethe ansatz and twisted transfer matrix 

To derive the master equation for the dynamical correlation functions we use the twisted transfer 
matrix in the framework of algebraic Bethe ansatz. The definition of this operator and related 
objects together with a brief sketch of their properties is given below. We refer the reader for 
more details to p. 

2.1 General framework 

The central object of algebraic Bethe ansatz is the quantum monodromy matrix. In the case of 
the XX Z chain Hl.lj) . this is a 2 x 2 matrix 




(2.1) 
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with operator-valued entries A,B,C and D which depend on a complex parameter A and act 
in the quantum space of states Ti of the chain. These operators satisfy a set of quadratic 
commutation relations given by the i?-matrix of the model. 

In the framework of algebraic Bethe ansatz, a sector 7-^(W2-a^) Qf ^j-^g space of states 7i with 
a fixed number N of spin down, = 0, 1, . . . , M, is spanned by vectors of the form 

N 

|V^) = n^(A,)|0), (2.2) 

where |0) is the state with all spins up and Ai, . . . ,Xm are arbitrary complex numbers. The 
dual space is constructed in terms of 

N 

(^| = (0|nC(A,), (2.3) 

i=i 

where (0| = |0)+. 

The Hamiltonian H^^^ of the homogeneous XX Z chain is given by the 'trace identity' (|1.5() 
in terms of the transfer matrix T(A) = ^(A) + D{X). All these transfer matrices commute, 
[T{\),T{^)] = 0, and their common eigenstates (which coincide with those of the Hamiltonian) 
are usually constructed in the form (|2.2j) (or (|2.3j) ) for Ai,...,Aiv satisfying a system of 
algebraic equations first obtained by Bethe in p . 

The local spin operators can also be expressed in terms of the entries of the monodromy 
matrix, by solving the quantum inverse scattering problem I.SOj : 

af = ■ tr(T(r//2)a-) ■ (2.4) 

Here cr° in the r.h.s. acts in the auxiliary space of T{X), while a" in the l.h.s. acts in the local 
quantum space Tij. 

2.2 Twisted transfer matrix 

It is actually convenient to introduce a slightly more general object: 
Definition 2.1. The operator 

T,(/x) = A{f,) + KDifx), (2.5) 

where k is a complex parameter, is called the twisted transfer matrix. The particular case of 
7^(/i) at K = 1 corresponds to the usual transfer matrix T{^). 

Let us consider the action of the operator Xt(^) in the subspace 7-^(^/2-^^) with fixed (but 
arbitrary) number of spins down A^. The eigenstates of Xi(/u) and their dual states in this 
subspace are denoted as |?/'k({A})) (respectively (•(/'^({A})!). They can be written in the form 
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()2.2j) (respectively (|2.3|l ). where the parameters Ai, . . . , A^v satisfy the system of twisted Bethe 
equations 

3^«(A,I{A}) = 0, j = l,...,N, (2.6) 

where 

N N 

ynifJ-lW) = aifi) Y[ sinh(Afc - fi + t]) + k d{fj.) sinh(Afc - n-rf), (2.7) 

k=l k=l 

and 

a(/i) = sinh*^(^ + r//2), d(/i) = sinh*^(^ - r?/2). (2.8) 
The corresponding eigenvalue of Xi(Ai) on IV'kII-^})) (or on a dual eigenstate) is 

( \{\\\ I \ A sinh(Afc - ^ + ??) . 1^ sinh(^ - A^ + ??) 

{A} = a(/i I I r-TTT ^ + M — 2.9) 

smh(Afc - \i) ^J-^ smh(// - Afc) 

Note that 

N 

3^«(/"l{A}) = n - • ^-(^l{A})- (2-10) 

fc=i 

Definition 2.2. ^ solution {A} o/ i/ie system h2. 6|) is called admissible if 

N 

d{\j) J] sinh(A, - Afc + ??) / 0, i = 1, . . . , A^, (2.11) 

and unadmissible otherwise. A solution is called ofF-diagonal if the parameters Ai, . . . , Aat are 
pairwise distinct, and diagonal otherwise. 

It is proven in Appendix A (see Theorem IA.1|) that the eigenstates corresponding to the 
admissible off-diagonal solutions of the system (|2.6I) form a basis in ■}i{m/2~n) ^ least if k is 
in a punctured vicinity of the origin (i.e. < |«;| < kq). 

Notations. Recall that in the particular case k = 1 the corresponding twisted transfer 
matrix is simply denoted by T in which the subscript k has been omitted. We follow the same 
agreement for all related objects. Namely, {ip,T,y) = ('(/'«; ''"kj 3^k) U=i- 

At K = 1, it follows from the trace identity (|1.5|) that the eigenstates of the transfer matrix T 
coincide with the ones of the Hamiltonian Hl.l|) . The corresponding eigenvalues can be obtained 
from (dig), (OH : 

N 



i/(0)|^({A})) = (X^ii;(A,)) • |^({A})), (2.12) 



where E{X) is called the bare one-particle energy and is equal to 



^^^^ - sinh(A+2Tsinh(A-2)- ^^-^^^ 
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One can similarly define the bare one-particle momentum. It is given by 

/sinh(A-2)\ 

Remark. The equation H2.12() holds if the parameters {A} of the eigenstate |'0({A})) corre- 
spond to an admissible off-diagonal solution of the system ()2.6|) at k = 1: 

N N 

3^(A, |{A}) = a{Xj) Yl sinh(Afc - A, + r?) + d{Xj) [J sinh(Afc - Xj - ??) = 0, j = l,...,N. 

k=l k=l 

(2.15) 

Unlike for generic k, the eigenstates basis of T (i.e. the eigenstates basis of the Hamiltonian) 
includes also states corresponding to unadmissible solutions of the system (|2.15() . We do not 
consider the expectation values of the operators with respect to such states. 

2.3 Scalar products 

We recall here the expressions for the scalar product of an eigenstate of the twisted transfer 
matrix with any arbitrary state of the form (|2.2|) or (|2..3j) . 

Let us first introduce some notations. We define, for arbitrary positive integers n, n' [n < n') 
and sets of variables Ai, . . . , A„, /ii, . . . , /i^ and vi, . . . ,Un' such that {A} C {v}, the nxn matrix 
^^iWdl^MM) as 



(^^K)jfc({A}, {/i}|{z^}) = a{^k)t{Xj,^ik) \\smh.{va - + ri) 

a=l 

n' 

- Kd{nk)t{fikAj) JJsinh(i/„ - ;Ufc - r/), (2.16) 
a=l 

with 

sinhri 

*A'^ = — T7^ T^-TTT r- 2.17 

^ smh(A - /i) smh(A - /x + ??) 

We also define, for arbitrary sets of variables Ai, . . . , A„ and ^i, . . . , ;U„, the Cauchy determinant 
XniilJ-}, {A}) as 

n 

n sinh(Aa - Xb) sinh(^fc - fia) 

-^n(M,{A})^det(-— Y^)=— n • (2-18) 

n Vsmh(/.,-A,); ^ sinh(;.,-A,) 

a,b=l 
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Proposition 2.1. [T^ Let Ai, . . . , Xjy satisfy the system \2. /xi, . . . , /xat 6e generic com- 
plex numbers. Then 

N N 

(oin^(/^^-)iV'K({A})) = (^«({A})in5(/^.)io) 

= n diXa) ■ X^'m, {A}) • det f _r,(;.,|{A}) j (2.19) 

a=l \ J / 

^_ diXa) detJ].({A},M|{A}). (2.20) 

n sinh(Aa - Afc) smh(/if, - fia) 

a>b 

The eigenstate |V'k({^})) is orthogonal to the dual eigenstate (^'^({/i})!, if the sets {A} and 
{/i} are different: {A} ^ {^} (see Appendix^. Otherwise 

(V^.({A})|^.({A})) = ^^^=^'^^^"^ • det {A}|{A}) (2.21) 

n sinh(Aa - Afc) 

a, 6=1 

= i-ir nli^Aj ^ / ^3;^(^^.||^})^ . (2.22) 



n sinh(Aa - Ab) 

a, 6=1 

The equations (|2.19j) - (|2.22|) are valid for any arbitrary complex parameter k, in particular 
at K = 1. In this case we denote = 0,^1^=1- 

3 Master equation for dynamical correlation functions 

Our main goal is to obtain an explicit representation for the time dependent correlation function 
of the third components of spin, 



K(0)<T^+l(t)) = ;'./r\.\i":.Tr^'r • (3-1) 



In this expression, |V'({-^})) denotes the ground state of the Hamiltonian (|l.ll) in the subspace 
j^{M/2-N) ^ Note however that all the results we present here for the finite chain remain valid 
for {A} being more generally any admissible off-diagonal solution of the system (|2.15p . 

Like in it is convenient to derive this correlation function from a special generating 
function. Indeed, let us consider the following time-dependent operator: 

Qt+^,m{t) = g) • (f ) • e^*^^°^ . g) . e-*^'"' , (3.2) 
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with Hji^ defined similarly as in (|1.5j) : 

2M cosh 77. (3.3) 



i?(°)=2sinh,^^7;-HA) 



^-2 



Using the expression (|2.4p for the reconstruction of the local spin operators in terms of the 
entries of the monodromy matrix and the fact that the twisted transfer matrix 7^ commutes 
with it is easy to see that 

xd(|).T— ^(|).e— , (3.4) 

_ 1 92 ^ 

= 2 ^^'^^'+1'™+! ~ Ql^+l,m ~ Q?+2,m+l + Q?+2,m)('^)\n=l- (3.5) 

Due to the translational invariance of the correlation functions, we can set I = and simply 
consider the following expectation value: 



In terms of this generating function, the time-dependent correlation function (|3.H) is thus given 

(3.7) 



K=l 



as ^ 

K(OX+i(t)) = 2(af (0)) -1 + 2DI ^Q,{m,t) 

where denotes the second lattice derivative defined as in (|3.5|1 . 

Like in the time- independent case, it is possible to derive a master equation for the generating 
function H3.6() in the finite chain. Indeed, we have the following result: 

Theorem 3.1. Let {Xi, . . . , X^} be an admissible off-diagonal solution of the system h2.15\) . 
Then there exists kq > such that, for < |k| < kq, the generating function QK,{m,t) 4,9. 6\) in 
the finite XX Z chain il.l]) is given by the multiple contour integral 

N r N 



N\ J 27ri 



r{±f }ur{A} 



X n^mh^(Aa-..)- ^ ^ r (3.8) 



n3^.(^a|{4)-det^(^^%P 



a=l 



/n i/lis expression, E(X) andp(X) denote respectively the bare one-particle energy and momentum 
\2.1'J\) and \2.14\l ; the integration contour is such that the only singularities of the integrand 
t^) within r{ib^} U r{A} which contribute to the integral are the points {±^} and {A}. 
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Remark 1. The master equation (|.S.8|) gives the expectation value Quirriji) with respect 
to an arbitrary eigenstate IV'II-^})) of ^ corresponding to any admissible off-diagonal solution 
{Ai, . . . , Xn} of (|2.15|) . In particular one can choose {A} such that IV'II-^})) is the ground state 
of the XX Z Hamiltonian. 

Remark 2. The master equation (|3.8() provides an integral representation of Qn{m^t) which 
is valid at least in a punctured vicinity of k = 0. On the other hand, it is clear from H3.2() and 
()3.6() that Qn^rriyt) is an analytical function of k for k £ C \ {0, oo}. Hence, the representation 
()3.8() can be analytically continued from any vicinity of the origin to the whole complex plane 
except K = 0, cxD. This does not mean, however, that one can set k to be an arbitrary specific 
value directly in the integrand of (|3.8j) . 

Remark 3. The time and space dependencies in the representation (|.S.8j) appear in a very 
suggestive way, namely through the exponent of the bare energy and momentum. Note however 
that they do not correspond a priori to any eigenstate since here Zj's are integration variables; 
nevertheless, the contribution of each point in the integration domain is measured in particular 
by the difference between the bare energy and momentum corresponding to the integration 
variables zj and the one corresponding to the ground state parameters Xj. 

As we already discussed in the Introduction, there are two possible ways to prove Theorem 

EH 

The first one, which follows the main strategy of PP, consists in acting with Ql^it) on 
the state (^/'({A})!, in computing the resulting scalar products and in rewriting the sum over 
partitions that follows as a single multiple integral of Cauchy type. In the next section, we 
explain how this approach, elaborated in can be applied in the time-dependent case in order 
to prove Theorem 13.11 

The second one was already announced in the conclusion of the master equation (|3.8|) 
can be obtained directly via a form factor expansion. Details are presented in Sectional 

The reader who is not interested in the technical details of these derivations can skip the 
next two sections and go directly to Section |H1 where a multiple integral representation of (|3.6() 
in the thermodynamic limit is deduced from p.8|l . 

4 Master equation via multiple action of transfer matrices 

In this section, we prove Theorem 13. II using the method developed in ^3] and pp. In order to 
apply this method, we need first to reduce the computation of the generating function 1)3. 6|) to 
the evaluation of the expectation value of some product of twisted transfer matrices. The idea 
is to reconstruct the operator ex.'p{itH^^) in 1)3.2(1 as a limit similar to (|1.6() : 

gii(//(°)+2Afcoshr,) ^ lim (TJ'l + e^-T-U^X)^ e = i2itsinhr?, (4.1) 

L—fOo \\2 / \1 J J L 
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and to use the fact that, for the specific value A = rj/2, the inverse operator^ is 
proportional to Xt(— r//2): 

T-\r^/2) = (Ka(r?/2)d(-r?/2))"V.(-r?/2). (4.2) 
This enables us to express the operator Qlmi^) ^ the limit 

«,„(t) = 2im Tig + .) .rr-(l) r'-d) .r-^(| + .), (4.3) 

Acting with the product of (untwisted) inverse transfer matrices on the eigenstate |'i/'({^}))) we 
obtain a representation of the generating function 1)3.6^ in terms of an expectation value of a 
product of twisted transfer matrices: 

Q«(m,t) = hm ^(-r?/2|{A})r-^(r?/2 + e|{A})(r,^(r?/2 + e)7;^— (-77/2)). (4.5) 

L— >oo 

It is convenient at this stage to introduce some arbitrary parameters uj\ , . . . , ^2L—m 

and to 

define (xi, . . . , X2L-m) = {uji + e, . . . ,ujl + e, wl+i - r/, . . . , u;2L-m - ??)• It follows that 



2L—m 2L—in 

Q.{m,t)= hm hm k'"-^ TT r-i(xa|{A}) ( TT T^)- (4-6) 



a=l a=l 



As soon as we have the representation 1)4. 6() . we can use directly the equation (4.15) of 
in order to obtain a multiple integral representation for the expectation value (O'^l^a))- 
leads to 



Q«(m,t)= hm hm ^'"-^^ / f] ^ 'jf T^pMl 
^ L^<^^.^1 N\ J ^}2m J-i r(rc,|{A}) 



r{x}ur{A} 

X TT /,r .^ •det^].(M,{A}|{.}) • f !'!'™ , (4.7) 
l}^yK{za\{z}) N det^J]({A},{A}|{A})' ^ ' 

where the closed contour T{x} U r{A} surrounds the points xi, . . . , X2L-m and Ai, . . . , Xn and 
does not contain any other singularities of the integrand^ . Formally the integrand in the equation 
()4.7p coincides exactly with its time-independent analogue in The difference, however, is 
that the number of parameters Xj eventually becomes infinite, and that one has to take the 

^One can easily obtain the equation l|4.2^ by applying the product 7^(?7/2)Xt(— 77/2) to an arbitrary state 

^More precisely, r{x} U r{A} is the boundary of a set of polydisks 'D{xa,r) and V{Xt,r) in C^. Namely, 
r{a;} — vJ^J^"^'V{xa,r), where 'D{xa,r) — {z € : \zk — Xa\ ~ r, k = 1,...,N}. The integration contour 
r{A} is defined in a similar manner. The radius r is supposed to be small enough. 
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homogeneous limit ujj — > r]/2 before the hmit L ^ oo. In this limit, some of the Xj go to 
T7/2 + e, while the remaining ones tend to —r}/2. Thus, to be able to proceed to this limit, one 
has to verify that all the solutions of the system 

N N 
yii{zj\{z}) = a{zj) sinh(2:fc - Zj + r]) + Kd{zj) sinh(zfc - Zj - rj) = (4.8) 

k=l k=l 

which contribute to the integral (|4.7|) are actually separated from the points r//2 + e, —rj/2 and 
the parameters {A}. Note that, thanks to the analyticity of {Qim) k € C \ {0,oo}, it is 
sufficient to evaluate this generating function in the punctured vicinity of k = 0. Observe also 
that one can set r//2 + e to be as close to 77/2 as necessary, since eventually we have to proceed 
to the limit L — > 00. 

The analysis of the solutions of the system H4.8|) performed in the article T was based on 
the results of (HH. However, most of those results were formulated in the case of inhomogeneous 
Bethe equations 

M N M N 

J|sinh*^(zj-^a + ??) sinh(2;fc-Zj+r/)-K sinh^(zj -^a) sinh(zfc - - ry) = 0, (4.9) 

a=l fc=i a=l fc=i 

for generic^ parameters ^i,...,^a/. We assume that all solutions Zj{K) of the homogeneous 
system i)4.8|) can be obtained from the solutions Zj{K\{^}) of the inhomogeneous system 1^4. y\) 
in the limit — > ry/2. Then, just like in the time-independent case, unadmissible and diagonal 
solutions of the system (j4.8|) do not contribute to the integral (|4.7j) , for det or det Q vanishes 
at these points, and the only solutions that actually contribute are admissible off-diagonal ones. 
Due to Lemma IA.H all of them are separated from the points ±^//2 (and hence from r]/2 + e as 
well, for £ small enough), at least in a vicinity of k = as far as k 7^ 0. For small enough, 
they are also obviously separated from the parameters {A}, which correspond to an admissible 
off-diagonal solution at k = 1. Thus, we can formulate 

Lemma 4.1. Let {A} he an admissible off-diagonal solution to \2.1f^) . and cjj — > 77/2. There 
exists kq > such that, forO < \k\ < kq, one can define a closed contour r{^}ur{—^}UT{X} 
which satisfies the following properties: 

1) it surrounds the points r//2, — r//2 and {A}, while all admissible off-diagonal solutions of 
the system \4-^ are outside of this contour; 

2) for L large enough, the only poles which are inside and provide non-vanishing contribution 
to the integral \4. 7| ) are Zj = r]/2 + e, Zj = —rj/2 and Zj = A^; 

3) for L large enough, the only poles which are outside (within a set of strips of width iir) 
and provide non-vanishing contribution to the integrand of \4- are the admissible off-diagonal 
solutions of the system \4-^ - 

■^Or at least 'well separated' (see |34) for definition). 
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Thanks to this Lemma, we can now take the hmit lvj ri/2 in (|4.7j) . It is easy to see that 



1- f-™'''rf ^T^l 

11 T{Xa\{\\) ) 



N 

n 



b=l 



sinh(z;, + ^ — e) sinh(Ab — 2 ~ ^)\^ ( sinh(z(, — ^) sinh(A;, + ^) 



sinh(zb 



ri^ \ L—m 

2' 



e) sinh(Ab + 3 — e) y \sinh(zb + 2) sinh(Af, — S) 



1 + 



sinh*^ e 



N 

n 



sinh(Afe 



3r, 



sinh^^(r7 + e) sinh(A6 + ^ - e) 



1 + 



K sinh e 



M ^ N 



n 



sinh(zb 



3»? 



e) 



sinh^(r? + e) ^j^^ sinh(z6 + ^ - e) 



(4.10) 



As the contour defined in Lemma [4.1l is independent of L, one can also safely take the limit L — > 
00 in the integrand. Using the definition of the bare one-particle energy E{\) and momentum 
p{\) (EUni), (EUl), one obtains 



hm hm ^--^'rf ^^T^Ml 



\ r{xa\{\})) 



N 



n exp (it{E{z^) - E{Xb)) + im{p{zb) - p{Xb))), (4.11) 



b=l 



which leads to the following multiple integral representation for the generating function Qn{m, t) 
in the finite chain: 



r{±^}ur{A} 



N 

n 



det^i.(M,{A}iM)- f!"^!!^M:!l!'!! - (4.12) 



fi3^«(^a|{4) ^ -vLJ'LJiLj; detjvf^({A},{A}|{A}) 
It remains to use p.l9|) - (|2.22|l and we obtain the master equation (|3.8|) . 



5 Master equation via form factor expansion 

The master equation (|3.8)) can be obtained via the form factor expansion for the expectation 
value ©^(m, t). This way is of course much shorter, for it is in particular not necessary to 
construct the generating function as a special limit of a product of twisted transfer matrices. 
However, the summation over form factors requires the completeness of the set of eigenstates 
of (see Theorem lA.ip . 

Inserting in 1)3. 6() the complete set of the eigenstates |'i/'K({/^})) of the twisted transfer ma- 
trices between the operators T^{^) • exp{itHj^^) and T~'"(2) • exp{—itH^^^) (see (|3.2() ). one 
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obtains 

where the sum is taken with respect to all admissible off-diagonal solutions {/ii, . . . , (J-n} (such 
that —tt/2 < 9(/ij) < 7r/2) of the system (|2.f)|) . As usual, we suppose small enough, but not 
zero. We have also used that 



i^p(Ab) = logT,(r?/2|{A}) -Mlogsinh?7, (5.2) 

6=1 

^^V.({A})) = • IMim, (5.3) 



which follows from ()2.9p . H2.14() and The scalar products in ()5.1() can be written in terms 

of Jacobians via ^^Tm - f^Tni : 



N N 



Q,{m,t) = {-1)^Y1 n sinh2(A„-/X6)-ne^*(^('^'')-^("^))+^"(^(^'')-pM 

a,b=l 6=1 



This last sum can be presented as a contour integral in 



n sinh^(A. - Z6) • ^ ^LL^^^ (5.5) 



dX 



-'"=1 n y^{Za\{z}) ■ detN ' MMM) 



a=l ^ ^ 



where the contour r{/x} surrounds all admissible off-diagonal solutions of the system ()4.8p and 
does not contain the points ±??/2 and {A}. The existence of such a contour follows from Lemma 
14.11 The factor 1/A^! appears due to the invariance of the off-diagonal solutions with respect to 
the permutations of {z}. 

Let us observe finally that the integrand in (|5.5j) is m-periodic with respect to each z^, and 
that it vanishes as soon as any z^ ±oo. Hence, the sum of the residues of the integrand 
within a set of strips of width iir is equal to zero. One can therefore evaluate the integral over 
the contour T{fi} by taking the residues outside this contour within this set of strips, i.e. inside 
the contour Fjit^} U r{A}. This leads to the master equation 1)3. 8() . 
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6 Thermodynamic limit 



Starting from this section the parameters Ai, . . . , Xn correspond to the ground state of 
the Hamiltonian in the subspace 7y(^^/2~^). This means in particular that 9Aj = if 9?r/ = 0, 
and 5ftAj = if Qr] = 0. Recall also that it is enough to consider the case h > 0, what implies 
N < M/2. 

In we showed how to obtain the thermodynamic limit of the generating function Qnim, 0) 
from the master equation by evaluating explicitly the integrals over r{A} and setting d{zk) = 
in the remaining integrals. A similar method applies to the time-dependent case, although 
the existence of the essential singularities at in the integrand makes this procedure more 

subtle. 

Let us start with the representation (|4.12|) . Due to the symmetry of the integrand with 
respect to the set {z} we have 

„ N N „ n r ^ 

f n^^i=E^^ / n^^^- / u (6-1) 

r{±|}ur{A} r{±f } ■'=^ r{A} i="+i 

In this expression, the integral over the contour r{A} can be rewritten as the sum over partitions 
of the set {A} into two disjoint subsets {A} = {A^^} U {Xa_}, with #{Aq,^} = n: 

N ^ „ n , n n 

1 X - / -i-r azi 



itE{Xi,)—imp(Xf,) 



a, =n 



detnn^{{z},{Xa^}\{z}u{Xa_}) detjvl7({A},{z}U{A»_}|{A}) 
U:=iy^i^a\{z}U{X^_}) ■ det^l7({A},{A,JU{A,_}|{A})' 

where the elements in the sets {z} U {Aa_} and {Aq^} U {Aq,_} are ordered accordingly. Here 
we have used that 



Res 

{z„ + l,...,Z]v} = {Aa_ } 



detni^{{zi, . . . ,zn}, {Xa+} U {Xa_}\{zi,. . .,zn}) 



N 

ny^iXalizi, ...,Zn}U {Xa_}) ■ dct , . . ■ , Zn} , {Xa,}\{zi , . . . , Zn} U {Aq,_}). (6.3) 

n 

Using now the system of Bethe equations (|2.15|) for variables {A} we can present the inte- 
grand of (|6.2j) as follows 



detnn,{{z},{K^}\{z}U{Xa_}) detA.l]({A},{z}U{A^_}|{A}) ^ , , . .... 
ULiy^i^a\{z}U{X^J) det^l7({A},{A,JU{A,_}I{A}) n '^^^ 

n 

^ b=l ag«+ detjv ^'{{z},{Xa_}\{X}) 

n sinh(A. - A, + r?) ft sinh(.. - + r?) det^ ^'({A« J U {A._ }) " 

a,fe6o+ a, 6=1 
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Here 



(Mt^)jk{{K+},U}) =t{Zk,\j)+Kt{\j,Zk) W 



and 



with 



'jk 



KiX) 



A=A, 



sinh(Aa - \j + 7?) -A- smh(Aj - Zg + r]) 
sinh(Aj — Xa + ri) ^J^ sinh(za - Xj + 77) ' 



N 



Y.K{Xj-Xa) 



a=l 



+ K{Xj-Xk), 



sinh 2r] 



(6.5) 



(6.6) 



(6.7) 



sinh(A — r]) smh(A + 7]) 

In $'({Aq,_|_} U {Aq,_}), the columns are ordered such that the n first ones correspond to the 
subset {Aq,_^} and the N—n last ones to {Aq,_ }. The matrix "if' has a more complicated structure. 
For A; > n its entries ^'^•^ coincide with the corresponding entries in the N — n last columns 
of $'({Aq, } U {Aq,_}). For the first n columns one has 



a{zk)t{Xj , Zk) 



N 



^1 



ci(..)t(z„A,)n -^iSir^^ 



a=l 



jk 



a{Zk_ 



n 



fe=i 



bee 



sinh(2fc-At,+?7) 
sinh(zfc-A(,-»7) 



k < n 



(6.8) 



In the thermodynamic limit M, ^ 00, M/N = const this expression can be simplified. 
Indeed, since d{z) and a{z) have zeros of order M at 2; = r]/2 and z = —r]/2 respectively, one 
can show that the contributions of the corresponding terms to the total result are bounded by 

/N\. Therefore, at M, N — s- 00, one can set d{z) = if 2; is in the vicinity of 77/2 and a{z) = 
if z is in the vicinity of — 7?/2. This gives us a simplified representation for the matrix elements 
(|6.S|) for A; < n in the thermodynamic limit: 



i(Aj, Zk)^ 



lim ^^-fc = 



M- 



jk 



Zk ^ ^ 



2' 



— K 



-1 



t{zk, Xj ^ 



n 

ago. 



sinh(zfc-Aa+}7) -pr sinh(zfc-Zft-?7) 
sinh(2j.-Aa-»7) 11 sinh(2;fc -25+77) ' 

0=1 



2 
2- 



In fact we can say that the limiting value ^'-^ has a cut between the points r//2 and —rj/2. 

The remaining steps are quite standard (see Jl], JH])- In the thermodynamic limit the 
distribution of the ground state parameters {A} can be described by the spectral density ptotW- 
In its turn the spectral density is a particular case of the 'inhomogeneous spectral density': 
PtotW = p{X, z)\^^,^^2- This inhomogeneous spectral density satisfies an integral equation 



-2mp{X,z) + / K{X 
Jc 



p)p{p,z)dii = t{X,z), 



(6.9) 



where the integration contour C in (|6.9() depends on the phase of the model. In the massless 
case —1 < A < 1 the contour C is an interval of the real axis [— A/j, A/j]. The boundary A/j 
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depends on the value of the magnetic field, in particular A/^ ^ oo at /i — > 0. For A > 1 (77 < 0) 
the limits ±A/i are purely imaginary, more precisely the integral in is taken over an interval 
of the imaginary axis. In particular A/^ = —iii/2 at /i = 0. 

In the thermodynamic limit, one can compute the ratio of the determinants det and det <!>' 
in terms of the inhomogeneous density. Indeed, since the last N — n columns of the matrix ^' 
coincide with the ones of the matrix ^'{{\a+} U {Aq,_}), we have 

yXi,Xa^ma+^) det^(<I>')^./*:J. (6.10) 



detTv $'({A„, }U{A„_}) 



.4 = 1 



Using the results of we have, for z in a vicinity of r//2. 



y i^')-h{h,z)^ ''/A , M^oo. (6.11) 



i=l 



If z is in a vicinity of —rj/2, we can set z = z — r] which gives t{z, Xj) = t{Xj, z), hence 



N 



V($0"'i(-,A.)^#^ = #^, M^oo. (6.12) 



Thus, in the thermodynamic limit M, N 00, N/M = const one has 



det^^;.,(M,{A„_}|{A„J) 



det[K(A„Zfc|{A„J,{z})], (6.13) 



detjv $jfc({Aa+}U{A«_}) 

where the function 7^^(A, 2:|{Ai, . . . , A„}, {zi, . . . , z„}) is defined differently in the vicinities of 
r]/2 and —ri/2: 

{p{X,z), 2~r?/2; 
- I n >^iMz~^b+v)>^iMzb-z+v) _ (6.14) 

-K p{A,Z + r]) [[ smh(Xt~z+rj)smh(z-z^+-n)^ Zr~.-7]/Z. 
6=1 

It remains to replace in (|6.2|) the sum over partitions of {A} by integrals over the support 
of the spectral density and we arrive at the multiple integral representation for the dynamical 
correlation function of the third components of spin: 



where 



(6.15) 

ft=i 



n (rr, -J— f I TT ^ TT sinh(Aa -Zb + 7]) sinh(zb - Ag + r?) 

t) - 2^ ^^^^ J ^ f ii 2m • \\ sinh(A, - A, + ^) sinh(z, - + r?) 

X TTei<(s(-5)-i?{A,))+im(p{-i>)-P(A,)) detM«({A},{z}) ■ det[K(A,-, Zfc|{A}, {z})]. (6.16) 
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The contour r{ib^} surrounds the points ±^ and does not contain any other singularities of 
the integrand. The parameter k in (|6.16j) is an arbitrary complex different from 0, oo. The 
functions entering the integrand are defined in (|2.1H|) . (|2.14|) . (|6.5|) . (|H.9j) and (|H.14|) . 

Due to the factors ex.Y>{itE{zb)) the integrand in 1)6. 16(1 has essential singularities in the 
points ib^. However, in the case i = 0, these essential singularities disappear and the integrals 
around — ^ vanish. The remaining part of the integrand has poles of order m at zj = ^. Hence, 
at t = the sum over n in ()6.16|) is actually restricted to n < m, and we reproduce the result 
of ^3] for the equal-time correlation function of the third components of spin. 

In Appendix El we explain how one can deduce from this expression the result obtained in 
|29j in the case of free fermions. 

Using exactly the same method one can obtain an integral representation for the corre- 
lation function for the partly inhomogeneous XX Z chain, where we associate a set of generic 
complex numbers ■ ■ ■ ,(,m, with the first m sites of the chain. In this case one should replace 
in the result ()6.16|) the one-particle momenta by their natural inhomogeneous modifications 

f\\ ^ V^i f sinh(A-Cfe) \ 

PinhW = — > log — — . . (6.17) 

^f^i Vsinh(A-Cfc + 7?)y 

The homogeneous case then corresponds to the limit = r]/2, k = 1, . . . ,m. 



7 Conclusion 

In this article, we have obtained a multiple integral representation for the dynamical correla- 
tion function. It is clear, however, that the method based on the master equation can be applied 
to other dynamical correlation functions as well. In fact we have seen that the time-dependent 
master equation for the generating function Qk(?7i, t) differs from its time- independent analogue 
only by the factors exp{it{E{z) — E{X)), which automatically appear in the framework of the 
form factor expansions for the correlation functions. It is quite natural to expect that the same 
simple modification holds also for other correlation functions. 

One interesting further development would be to obtain a generalization of the multiple 
integral representations for the dynamical correlation functions at finite temperature. A method 
to consider temperature correlation functions by algebraic Bethe ansatz was proposed recently 
in j2L)j . It is possible that this technique can be successfully combined with the approach used 
in this paper. In particular one obvious question is whether there exists also a master equation 
for the temperature-dependent case. It would lead also to the interesting question of the form 
factor expansion at non-zero temperature. 

It is also well known that, for the case of free fermions A = 0, the dynamical correlation 
functions of the XX Z chain satisfy difference-differential classical exactly solvable equations 
|35l 1361 07] . It is natural to wonder whether this property holds also for general A, or at least 
for some specific cases. We hope that the multiple integral representations for the dynamical 
correlation functions open a way to study this problem. 
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A Admissible solutions of the twisted Bethe equations 

Let e^^-' = Xj and = q. Then the system (|4.8|) takes the form 



N N 



Y^{xj\{x]) ^ {x, - q-^Y' \{{x, - q^Xa) - Kq^""-^-^' (x,- - qY'J{{x, - q-^Xa) = 0. (A.l) 



a=l a— 1 



It is clear that, in the limit k ^ 0, all admissible solutions of HA.1|) go to q~^ and the Jacobian 
matrix {dYf-{xj\{x})/dxk) has the rank zero at k = and Xj = q~^. Our goal, however, is to 
solve the system HA.1|) for \k\ small enough, but not zero. 

A simple example shows that, unlike in the inhomogeneous case, the solutions xj^k) of ()A.1|) 
are not holomorphic functions at k = 0. Indeed, let q = i (free fermions). Since for admissible 
solutions Xj + Xk ^ 0, we obtain 



Xj + I 



K e 



i^(N-i-M/2)^ for q = i. (A.2) 



This system has an obvious solution 

+ 1 . . ., , 1 ITT 



-, where 0^ = -i|K|Met(2"j+^-i), € {0, 1, . . . , M - 1}. (A.3) 



1 



Thus, in this {km ) and different choices of the branch of km correspond to different 

solutions. 

One can treat the general case similarly. 

Lemma A.l. There exists kq > such that, for < |k| < kq, all admissible solutions of the 
system \A . 1\) belong to the vicinity of q~^ , but are separated from this point. 

Proof. Let us make in ()A.1|) the substitution Xj = q^^ + 6uj, where 9 is one of the solutions 
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of 9^^ = Kq^^ ^ *^ . Then we arrive at 

N 

Ye{uj\{u]) = uf W (e{uj - q^ua) +q-q-' 



a = l 



N 



(9uj - q + q 



-1\M 



q "^Ua) + q ^ 



,-3 



0. (A.4) 



At = one has 

%(0) = {q-^-q) 

and 

dYe{uj\{u}) 



2(1-JV) 2iTvnj 
\ M . e M 



duk 



Uj G {0,1,...,M-1}, 
MS,k-uf-\0)-{q-q-Y''- 



(A.5) 



(A.6) 



9=0 



Hence, due to the imphcit function theorem in a vicinity of = 0, there exists a unique 
holomorphic solution Uj{9) of the system ()A.4p which takes the value ()A.5|) at = 0. This 
implies that 

xj{k) = q-^ + 9uj{0) + o{9). (A.7) 

Therefore, for \k\ small enough, but not zero, all the admissible solutions of the system ()A.1|) 
are in a vicinity of but Xj{K) q^^ ■ It is also evident that the replacement 9^9' = 9e m 
corresponds simply to a different choice of the integers nj in ()A.5|) . □ 

Observe that, due to HA.7|) . the admissible solutions of the system ()A.1|) for \k\ small enough 
are off-diagonal if the integers nj in (|A.5|) are pairwise distinct. Thus, the system HA.1|) (and 
therefore the system ()4.8() ) has = dim(7^(*^/-^~^)) different admissible off-diagonal solutions. 
In the following two lemmas we prove that the corresponding eigenstates \ipK,{{z})) are linearly 
independent. 

Lemma A. 2. Let the sets {z} and {z'} be two different admissible off-diagonal solutions of the 
system \4-^ for k ^ 0. Then 

= 0. (A.8) 

Proof. The scalar product (|A.8|) is given by (|2.2fl|) . If 7^ Zj, Vj, k, then one can express 
Kd{zk) in terms of a{zk) via (|4.8|) for the parameters {z}. Then we obtain 



N 



b=l 



N 

n sinh(4 - Zfe + ?]) 

a,b=l 

~N AT 
n smh(2;a - Zb) smh(2;^ - z'J 

a>b 



detM({4,{z'}), (A.9) 
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where M{{z}, {z'}) = Mfs({z}, {z'}) (see (|6.5|) ). It was proved in ^S] that this matrix has 

K = l 

an eigenvector with zero eigenvalue 



N n n 

^ MjkVk = 0, where ^^ = 11 sinh(4 - Za) JJ sinh~^(4 - 4)- (A.IO) 

k=l a=l a=l 

Thus, in this case, the scalar product vanishes. If several parameters z' coincide with z, say, 
Zj = Zj for j = 1, . . . , n, then one should first proceed to this limit in the first n columns of the 
matrix 0^ (|2.2fl|) . and only in a second step use the equations (|4.8|) . One can easily verify that 
in this case the matrix we obtain has a zero eigenvector of the same type as Vk- Q 

Lemma A. 3. There exists kq > such that, for < [k| < kq, 

{M{z})\M{z}))^o- (A.ii) 

Proof. The 'square of the norm' of \iPk{{z})) is proportional to the Jacobian (|2.22j) . The 
last one coincides with the Jacobian (|A.6|) up to a trivial factor which do not vanish at k ^ 0. 
Since dYg{uj\{u}) / duk is a continuous function of 9 and due to ()A.6|) . it is also non- vanishing 
in a vicinity of k = 0. □ 

Using these two lemmas we prove 

Theorem A.l. There exists kq > such that, for < |k| < kq, the states \iPk{{z})) corre- 
sponding to the admissible off-diagonal solutions of the system i4-^ form a basis in the subspace 

^(M/2-7V) _ 

Proof. It follows from Lemmas IA.21 IA.3I that 

1) Different admissible off-diagonal solutions of the system ()4.8I) correspond to different 
states \^pK{{z})). Hence, the total number of the last ones is dim(7^(*^/2~^)). 

2) These states are linearly independent. □ 

B Free fermions 

In the particular case of free fermions A = (i.e. rj = — i7r/2), the dynamical correlation 
function (af (0)o"^^_]^(t)) can be written in the form ^Hj : 

2 '^^ 

(a^(0)a^+l(t)) = - 1^ +^ j e^^*™^P+*™f dp j ^-i^t cos q+imq 

~kp [— 7r,7r]\[— fcj?,fcj?] 

Here cos kp = h/A. This result was obtained by means of a summation over the complete set 
of excited states. The constant term in ()B.1|) corresponds to the diagonal contribution (cr^)'^. 
The integration variable p belongs to the Fermi sphere [—kp, kp] and describes the ground state 
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distribution. The integral over the parameter q corresponds to the thermodynamic limit of the 
sum over the excitations outside the Fermi sphere. 

In this appendix, we show how to reproduce this result by the direct use of the representation 
(|™) . 

At r] = — z7r/2 the integral equation 1)6. 9(1 is explicitly solvable: the inhomogeneous spectral 
density has the form 



p{X,z) 



(B.2) 



vr sinh 2(A — z) ' 

whereas the boundary of integration Ah is defined by the relation cosh 2 Ah = 4:/h. The function 
TZ'^ at 7] = —in/ 2 becomes 



KiX,z\{X},{z}) 



z ~ 



-iTr/i; 



7rsmh2(A-2) ' 



(B.3) 



The main simplification, however, comes from the matrix M^, which at A = is proportional 
to «; — 1: 

2(k - 1) 



(M,U({A}|{z}) 



(B.4) 



sinh 2{Xj — z^) 

Thus, after taking the second derivative over k and setting k = 1, all the terms of the series 
(|6.16|) with n > 2 vanish. 

(2) 

Let us first consider the term [m,t) corresponding to n = 2. After differentiating with 
respect to k one has 



1 



327r4 



d^X 



cfz ■ det 



^itE(zk)+imp(zk) 

sinh(Aj — Zk) 



det 



r{Tf } 



.—itE{\j)—imp{\j ) 

sinh(Aj — Zf^) 



(B.5) 

Like in the time-independent case the integral over z^ can be taken by the residues outside of 
the contour r{=F^}, i.e. in the points Zk = Xj. This gives us 



K=l 



I cfXdet 



I _ gii{S(A,)-£;(Afc))+im(p(A,)-p(Afc)^ 

sinh(Aj — Afc) 



(B.6) 



It remains to take the second lattice derivative and, after the standard change of variables 
cosh2Aj = cos~^pj, we obtain 



2Dl-^Q'i\m,t) 



K=l 



Akl 1 



g4ii cos p+imp 



kp 



(B.7) 
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The term Ql^''(m,t) corresponding to n = 1 appears to be more complicated. We have 
Ah / 



+ K 



-1 



dz 



-Ah r{f} 
Evaluating the integral over r{— ^} we obtain 
dz 



sinh (A — z) 



Jt{E{z)-E{\))+im{p{z)-p(\)) 



(Bi 



r{^} 



sinh^(A — z) 



. it{E(z)-E{\))+im(p{z)-p(\)) 



dz 



r{f } 



sinh (A — z] 



. ^t{E(,z)-E{\))+rm(,p{z)-p{,X)) ^ 27r[t^'(A) + mp' [X)]. (B.9) 



It is clear that the second lattice derivative of the last term vanishes, and we have 

Ah 



Di / dX 



dz 



-Ah r{f} 



sinh (A — z) 



. ^it{E{z)-E{X))+im{p{z)-p{X)) _ (B.IO) 



We can now explicitly differentiate this expression with respect to k and m, which leads to 

Ah 



K=l 



/ dX 



dz 



cosh 2A cosh 2z 



. ^it{E{z)-E{X))+imip{z)-piX)) _ (B.ll) 



-Ah r{if} 

It remains to move the contour r{^} to the boundaries of the strip 



dz = / dz 



dz. 



(B.12) 



r{f } 



and after the same change of variables as in the case n = 2, we finally obtain 



d 



1 

«=i J-k 



^iit cos p+imp 



dp / e 



-iit cos q+imq 



dq. (B.13) 



Using now that (cr^) = 1 — 2kp/TT we reproduce the result ()B.l 



*In fact one can obtain (cr^) form the same generating function QK(m, i) by taking the first derivatives with 
respect to m and k at k = 1. 
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